Conductance enhancement in quantum point 
contact— semiconductor— superconductor devices 
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We present numerical calculations of the conductance of an interface between a phase-coherent 
two-dimensional electron gas and a superconductor with a quantum point contact in the normal 
region. Using a scattering matrix approach we reconsider the geometry of De Raedt, Michielsen, and 
Klapwijk [Phys. Rev. B, 50, 631 (1994)] which was studied within the time-dependent Bogoliubov- 
de Gennes formalism. We find that the factor-of-two enhancement of the conductance G N s compared 
to the normal state conductance G N for ideal interfaces may be suppressed for interfaces with a 
quantum point contact with only a few propagating modes. The suppression is found to depend 
strongly on the position of the Fermi level. We also study the suppression due to a barrier at the 
interface and find an anomalous behavior caused by quasiparticle interference. Finally, we consider 
the limit of sequential tunneling and find a suppression of the factor-of-two enhancement which 
may explain the absence of conductance enhancement in experiments on metal-superconductor 
structures. 
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I. INTRODUCTION 

Charge transport through a normal conductor- 
superconductor (NS) interface is accompanied by a con- 
version of quasiparticle current to a supercurrent. In the 
Andreev reflection, by which the conversion occurs, an 
electron-like quasiparticle in the normal conductor (with 
an excitation energy lower than the energy gap of the 
superconductor) incident on the NS interface is retrore- 
flected into a hole-like quasiparticle (with reversal of its 
momentum and its energy relative to the Fermi level) and 
a Cooper pair is added to the condensate of the super- 
conductor For an ideal NS interface, the signature 
of Cooper pair transport and the Andreev scattering is 
a doubling of the conductance compared to the normal 
state conductance. 

A theoretical framework for studies of the scatter- 
ing at NS interfaces is provided by the Bogoliubov-de 
Gennes (BdG) formalism 0] where the scattering states 
are eigenfunctions of the BdG equation 



semiconductor heterostructure and a superconductor ex- 
perimentally, has provided a playground for investigat- 
ing the interplay between the Andreev reflection and 
the mesoscopic effects seen in mesoscopic semiconduc- 
tor structures In the recent years the technological 
efforts have revealed a variety of new mesoscopic phenom- 
ena, see Refs. and references therein. One class of 
the studied devices are the quantum point contact (QPC) 
2DEG-S and S-2DEG-S devices with the QPC in the nor- 
mal region. The dc Josephson effect and the quantiza- 
tion of the critical current in QPC S-2DEG-S junctions 
have been studied extensively both experimentally by e.g. 
Takayanagi and co-workers [|3|_ and theoretically by e.g. 
Beenakker and van Houten 
Furusaki, Takayanagi, and TsuE 
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which is a Schrodingcr-like equation in electron-hole 
space (Nambu space). Here 7i(r) is the single-particle 
Hamiltonian, A(r) is the pairing potential of the super- 
conductor, E is the excitation energy, and u(r) and v(r) 
are wave functions of electron-like and hole-like quasipar- 
ticles. 

The technological possibility of studying the interface 
between a two-dimensional electron gas (2DEG) in a 



The linear-response conducting properties of QPC 
2DEG-S structures have been studied by several groups. 
In the analytical work of Beenakker [|l3| a ballistic nor- 
mal region with a QPC modeled by a saddle-point po- 
tential was considered. The effect of elastic impurity 
scattering was considered numerically by Takagaki and 
Takayanagi [Q who considered a disordered region be- 
tween a narrow-wide (NW) constriction and the super- 
conductor. Both of these studies of the conductance were 
based on a scattering matrix (S-matrix) approach and 
the BdG formalism. In the numerical simulations of De 
Raedt, Michielsen, and Klapwijk |l5), based on the time- 
dependent BdG equation, a wide-narrow-wide (WNW) 
constriction was considered. Here, the aim was to study 
the electron-hole conversion efficiency and the robustness 
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of the back-focusing phenomena of the Andreev reflec- 
tion. 

One of the properties of the QPC is that most trans- 
mission eigenvalues are either close to zero or unity. For 
an ideal QPC 2DEG-S interface, the Andreev reflection 
will therefore give rise to a factor-of-two enhancement 
of the conductance G NS compared to the normal state 
conductance G N jl0| , |l3]| , which is quantized in units of 
2e 2 /h However, as pointed out by van Houten and 
Beenakkcr |l(J, deviations from the simple factor-of-two 
enhancement should be expected when the position of 
the Fermi level does not correspond to a conductance 
plateau. The presence of impurity scattering in the nor- 
mal region and/or interface roughness will also suppress 
the doubling of the conductance H] . 
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In the sequential tunneling limit the conductance can 
be found by considering the QPC and the interface as 
two series-connected resistive regions and in the limit 
W 3> W the enhancement of the conductance compared 
to the normal state conductance vanishes even for ideal 
NS interfaces. This may be an explanation for the un- 
expectedly low conductance enhancement in the experi- 
mental results of Benistant et al. jl8| on Ag-Pb interfaces 
where the current is injected into the Ag crystal through 
a point contact. 

The text is organized as follows: In Section II the S- 
matrix formalism is introduced, in Section III we formu- 
late our model, in Section IV the scattering scheme of 
the considered geometry is presented, and in Sections V 
we present results of several applications of our scattering 
scheme. Finally, in Section VI discussion and conclusions 
are given. 
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FIG. 1. Geometry of a WNW 2DEG-S junction with a 
hard-wall confining potential and a barrier at the 2DEG-S 
interface. 

Using an S-matrix approach, we study the linear- 
response regime of a phase-coherent ballistic QPC 2DEG- 
S system where the QPC is modeled by a WNW con- 
strictions with a hard-wall confining potential, see Fig. 
|l|. We report new results for the device studied by De 
Raedt et al. |n| which had a relative width W/W = 
1.7^m/(1.6 x 335 A) ~ 31.72, an aspect ratio Li/W = 
5/1.6, and a relative length L 2 /W = 20/1.6. Apply- 
ing the S-matrix formalism instead of the computation- 
ally more complicated time-dependent BdG formalism, 
we are able to study a larger part of the parameter-space 
where we also consider a barrier (with a normalized bar- 
rier strength Z) at the NS interface. We focus on the 
regime with only a few propagating modes in the QPC. 
In this regime the transmission eigenvalues of the QPC 
depend strongly on the actual position of the Fermi level. 
Even for an ideal interface this gives rise to a strong sup- 
pression of the conductance for certain positions of the 
Fermi level as predicted by van Houten and Beenakkcr 
p0{ and subsequently seen in the work of Beenakker |l^, 
Fig. 1]. In the presence of a barrier at the interface, the 
QPC gives rise to an enhanced tunneling through the 
barrier (compared to the case without a QPC) as in the 
case of rcflectionless tunneling effect of diffusive junctions 

@0- 



II. SCATTERING MATRIX FORMALISM 

The scattering approach to coherent dc transport in su- 
perconducting hybrids follows closely the scattering the- 
ory developed for non-superconducting mesoscopic struc- 
tures, see e.g. the text-book by Datta fl<| . 

For an ideal NS interface, the interface acts as a phase- 
conjugating mirror within the Andreev approximation jlj] 
and the rigid boundary condition for the pairing poten- 
tial 



A(r) = A e^e(x - L) 



(2) 



where A is the BCS energy gap |20|| , ip is the phase of 
the pairing potential, Q(x) is a Heaviside function, and 
L = L\ + Li is the length of the normal region (see Fig. 

0>- 

In the linear-response regime in zero magnetic field, 
Beenakker jl]| found that the conductance G = dl/dV 
is given by 
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(4) 



is a non-linear function of the transmission eigenvalues 
T n (n = 1,2, ... ,N) of itt. Here t is the N x N trans- 
mission matrix of the normal region, N being the number 
of propagating modes. 

Equation (|3|) holds for an arbitrary disorder potential 
and is a multi-channel generalization of a conductance 
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formula first obtained by Blonder, Tinkham and Klap- 
wijk |p2j who considered a delta function potential as a 
model for the interface barrier potential. The computa- 
tional advantage of Eq. (||) over the time-dependent BdG 
approach of De Raedt et al. M| is that we only need 
to consider the time-independent Schrodinger equation 
with a potential which describes the disorder in the nor- 
mal region, so that we can use the techniques developed 
for quantum transport in normal conducting mesoscopic 
structures. 

For more details on e.g. finite bias and/or tempera- 
ture, see Lesovik, Fauchere and Blatter p3| , Lesovik and 
Blatter , and the reviews of Beenakker || , and Lam- 
bert and Raimondi [ffl|. 



relevant for the geometry shown in Fig. |l|. The S-matrix 
S of a scattering region relates the incident current am- 
plitudes af to the outgoing current amplitudes bf. For 
a scattering region with two leads, S is a 2 x 2 block- 
matrix with sub-matrices Sn, S12, S21, and 522, where 
the diagonal and off-diagonal sub-matrices are reflection 
and transmission matrices, respectively. The appropriate 
scattering scheme for three scattering regions (the WN 
and the NW constrictions, and the interface barrier po- 
tential, respectively) connected by ballistic conductors is 
shown in Fig. |^. 

The WN constriction is described by the S-matrix 5 WN 



b+ \ ( o+ 
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III. MODEL 



We describe the geometry of Fig. [I] by the Hamiltonian 



W(r) = -— V + Vi(r)+V (r)-M ; 



(5) 



where fi is the chemical potential. The barrier potential 
is given a Dirac delta function potential with strength H 
P2|: 



the narrow region of length L\ by the propagation-matrix 

U N 



at)- U b+ 



the NW constriction by the S-matrix 



(10) 



(11) 



V,(t) =H6(x-L) 



(6) 



and the transverse motion is limited by a hard-wall con- 
fining potential 



, . / , \y\<W(x)/2 
c[ > \ 00 , \y\ > W(x)/2 ' 



(7) 



where the width W(x) defines the WNW constriction and 
is given by 



W , x<0 
W{x) = { W , < x < L x 
W , x > L x 



(8) 



The scattering states can be constructed as linear com- 
binations of the eigenstates of the Schrodinger equation. 



IV. SCATTERING SCHEME 

In the following subsections, we consider the S- 
matrices of a system with the Hamiltonian in Eq. (pT 



the wide region of length L2 by the propagation-matrix 



U v 



b± 



and the delta function barrier by the S-matrix S s 
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To apply Eqs. @ and (f§ we need to calculate the com- 
posite transmission matrix t = S21 which is a sub-matrix 
of the composite S-matrix S = S WN ®U N ®S NW ®U W ®>S S 
relating the outgoing current amplitudes to the incoming 
current amplitudes, 
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The meaning of the symbol ® is found by eliminating the 
internal current amplitudes |l9| . As a final result we find 
the transmission matrix 
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FIG. 2. Scattering scheme appropriate for the normal re- 



A. Quantum point contact 

We consider a QPC which we model by a WNW con- 
striction defined by a hard-wall confining potential, see 
Fig. [I]. This geometry has been considered by Szafer and 
Stone j25| and Weisshaar, Lary, Goodnick and Tripathi 
p6[ in the context of conductance quantization of the 
QPC in a 2DEG, and recently by Kassubek, Stafford and 
Grabert p7| in the context of conducting and mechani- 
cal properties of ideal two and three-dimensional metallic 
nanowires. We follow Kassubek et al. and calculate 
the composite S-matrix 5 WNW = 5 WN ® U N <g> S* NW . In 
zero magnetic field, where all S-matrices satisfy S = S T , 
the individual S-matrices are given by 
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g NW 
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(16) 
(17) 
(18) 



where X™ n , = 8 nn i exp (ik n Li) describes the narrow re- 
gion with free propagation of propagating modes and 
an exponential decay of evanescent modes. Here k n = 

fcp^/l — {niT /k-pW') 2 is the longitudinal wave vector of 
mode n in the narrow region. The S-matrices of the WN 
and NW constrictions are related through an exchange 
of leads. 



By elimination of the internal current amplitudes we 
find the composite transmission matrix 



' WNW L WNW 



^WNW ' WNW 



(19) 



where 



' WNW — ' NW 



l-(A N r NW ) Z X N r NW X N 4 w , (20) 



/ - / 1 - (X N r ) 2 } 1 X N t T 



(21) 



The S-matrix of the NW constriction can be found 
from a matching of scattering states which are eigen- 
states of the Schrodinger equation with the Hamiltonian 
in Eq. (||) where we only consider the part of the poten- 
tial which sets up the NW constriction. From a matching 
of scattering states we find that 



r»w = (1 + QQ ) 
t NW = 2g T (i + gg T ) 

r Nw = (e T Q + i) 1 f 

elements of the 
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(22) 
(23) 
(24) 



g-matrix can 
$ U1 ) where (4 



be 



writ- 



where the 

ten as g nw = ^/K w /k n {<f> n 
I-oc 4'n(y)^w(y) is an overlap between transverse 
wave functions of mode n in the narrow region and mode 

w in the wide region. Here K w — k-pyj 1 — (wir/rjk-pW') 2 
is the longitudinal wave vector of mode w in the wide 
region and r\ = W/W is the relative width of the con- 
striction. 

The overlap can be calculated analytically since its el- 
ements consist of overlaps between transverse wave func- 
tions 4> n and & w which are either two Sine or two Cosine 
functions (the overlap between a Sine function and a Co- 
sine function or vice versa is zero due to the odd and even 
character of the two functions). From the overlap-matrix 
we get the following elements of the g-matrix 
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(25) 



where the parity V(j) of j is V(j) 
V(j) ee -1 if j is odd. 



1 if j is even and 



In the numerical evaluation of Eqs. (22)-([24|) and Eqs. 
( pp| ) and (|l]) it is crucial to let the number of modes 
in the narrow and wide regions extend over both propa- 
gating modes and evanescent modes. After all matrix 
inversions are performed, the reflection and transmis- 



sion matrices are projected onto the propagating modes. 
In practice, numerical convergence of the reflection and 
transmission matrices is found for a finite cut-off in the 
number of evanescent modes. For the considered device, 
the number of evanescent modes is roughly ten times the 
number of propagating modes in the wide region corre- 
sponding to 1000 x 1000 matrices. 
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In the limit W 3> W, Szafer and Stone (||] employed 
a mean- field approximation for the overlap (<f> n \ <& w ) in 
which mode n of the narrow region couples uniformly to 



only a band of modes (with the same parity as mode n) 
in the wide region within one level spacing so that the 
elements of the g-matrix take the form 



/ ( k w W' 



1/4 



1-1/2 



(n — 1)?/ < w < (n + 1)?/ 
otherwise 



(26) 



Within this approximation, there is no mode mixing and 
tt^ becomes diagonal with the transmission eigenvalues 
along the diagonal. This approximation was found to 
capture the results of an exact numerical calculation ]25| ] 
when used with the Landauer formula, Eq. (|J), which 
is linear in the transmission eigenvalues T n . However, 
for an NS interface, the conductance formula, Eq. (||), 
is non-linear in T n which also makes off-diagonal compo- 
nents in ttf important. As we shall see (in Fig. |J), the 
mean-field approximation cannot reproduce the results 
of an exact numerical calculation of G NS in the same nice 
way as for G N . 



B. Wide region 

The wide region L\ < x < L is described similarly to 
the narrow region by 



1 + iZj COs9 v 

—iZ ' I cos 9 W 
1 + iZj cos# u 



(29) 
(30) 



where the normalized barrier strength is given by Z = 

H/hvp and cos9 w = K w /k? = yl — (wir / kpW) 2 . The 
results differ from those of a one-dimensional calcula- 
tion |22j since we have taken the parallel degree of 
freedom into account. However, if we introduce an 
angle dependent effective barrier strength Z etl (6 w ) = 
Zj cos 9 W ^,^] , the transmission and reflection ampli- 
tudes can formally be written in the one-dimensional 
form of Ref. [^2). The transmission eigenvalues of t s t\ 

in contrast to the 



are given by T£ = (l + Z% t (9 w )) 

mode-independent result = (l 
dimensional calculation 



of a one- 



x^ 







(27) 



where X™ w , = 5 WW > exp (iK w L 2 ) describes both the free 
propagating modes and the exponential decay of the 
evanescent modes. Here K w is the introduced longitu- 
dinal wave vector of mode w in the wide region. 



C. Interface barrier potential 

We consider an NS interface of width W with a bar- 
rier which we model by a Dirac delta function potential, 
following Blonder, Tinkham and Klapwijk p^ |. The S- 
matrix elements for the delta-function potential is found 
from a matching of scattering states which are eigenstates 
of the Schrodinger equation with the Hamiltonian in Eq. 
(|J) where we only consider the part of the potential con- 
sisting of the barrier at the interface. In zero magnetic 
field one finds the symmetric result 



S 5 = 



(28) 



with 



V. RESULTS 

A. Phase-coherent junction with ideal interface 

For the case of coherent transport through an ideal 
2DEG-S interface with a WNW constriction in the nor- 
mal region, see lower insert of Fig. ^, the conductance 
G NS and the normal state conductance G N can be found 
from Eqs. (|) and (|) with the transmission matrix 

t = twNW- 

Figure || shows the conductance as a function of 
k-pW /ir based on a numerical calculation of £ WNW (full 
lines) and the mean-field approximation (dashed line) for 
a WNW constriction with an aspect ratio L\jW' = 1 
and a relative width W/W = 31.72. The conductance 
G NS is seen to be approximately quantized in units of 
Ae 2 /h which is twice the unit of conductance for the 
normal state conductance G N . However, just above the 
thresholds (kpW' /ir = 1, 2, 3, . . .), oscillations due to res- 
onances in the narrow region of the constriction are ob- 
served. In the normal state result, these resonances are 
small but in contrast to the Landauer formula, G NS is not 
linear in the transmission eigenvalues and this makes the 
resonances much more pronounced compared to those in 
the normal state conductance. Another signature of the 
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non-linearity of G NS and the importance of off-diagonal 
transmission, is that the mean-field approximation is in 
good agreement with the numerical calculation for G N 
whereas it has difficulties in accounting for G NS . The 
sharpness of the resonances is to a certain extent due to 
the wide-narrow-wide constriction, and is suppressed in 
experiments with split-gate-dcfined constrictions. How- 
ever, as shown by the simulations of Maa0, Zozulcnko, 
and Hauge Q resonance effects do persist even for more 
smooth connections of the narrow region to the 2DEG 
reservoirs. 




o 



2 2.5 

kpW'/lT 

FIG. 3. Conductance G N s and normal state conductance 
G N of a coherent WNW 2DEG-S junction as a function of 
k^W /it. The constriction has an aspect ratio Li/W' = 1 and 
a relative width W/W' = 31.72. The full lines are for a numer- 
ical calculation of i WNW and the dotted lines are results within 
the mean-field approximation of Szafer and Stone J2f|. The 
upper insert shows the normalized conductance g = Gns/Gn 
also as function of kpW /n. 

The normalized conductance g = G NS /G N , shown in 
the upper insert, is two on the conductance plateaus 
but for certain "mode-fillings" of the constriction it is 
strongly suppressed and for kpW'/ir ~ 2 (two propagat- 
ing modes) we get g ~ 1.5. This effect, which occurs 
at the onset of new modes, was also seen in the calcula- 
tions of Beenakker (l^, Fig. I]. As the number of modes 
increases, these dips vanish and the normalized conduc- 
tance approaches its ideal value of two. The reason is sim- 
ple: suppose the constriction has N propagating modes, 
then the N — 1 of them will have a transmission of order 
unity and only a single mode (corresponding to the mode 
with the highest transverse energy) will have transmis- 
sion different from unity. As N increases, the effect of 
the single mode with transmission different from unity 
becomes negligible for the normalized conductance and 
from Eqs. (Jfy and (Q) it follows that liniTv-^oo.9 = 2. 

Since the quasiparticle propagation is coherent and the 
Andreev scattering is the only back-scattering mecha- 
nism, the phase conjugation between electron-like and 



hole-like quasiparticles makes the conductance G NS inde- 
pendent of the separation Li of the constriction and the 
interface. If evanescent modes in this region were also 
taken into account the results would depend weakly on 
L2 as it was found in the simulations of De Raedt et al. 
|l5f and our results should be compared with their results 
in the large L2-limit. As we shall see below, interfaces 
with a finite barrier (and thereby normal scattering at 
the interface) lead to size-quantization and thereby reso- 
nances which will depend on L 2 . 

The back-focusing phenomenon of the Andreev re- 
flected quasiparticles and the lowering of the normalized 
conductance due to a QPC in the normal region was stud- 
ied by De Raedt, Michielsen and Klapwijk |15j by solv- 
ing the time-dependent BdG equation fully numerically. 
In their wave propagation simulations, the QPC is also 
modeled by a WNW constriction with a relative width 
W/W = 1.7|Um/(1.6 x 335 A) ~ 31.72, an aspect ratio 
Li/W = 5/1.6 and a relative length L 2 /W = 20/1.6. 
For the particular "mode-filling" tpW/n — 3.2, they 
find a normalized conductance g = 1.87 < 2, but the de- 
pendence on the "mode-filling" was not studied in detail. 
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FIG. 4. Normalized conductance g = G N s/G N of a coher- 
ent WNW 2DEG-S junction as a function of kwW'/ir. The 
constriction has an aspect ratio Li/W' = 5/1.6 and a rela- 
tive width W/W' — 31.72. The curve is based on a numerical 
calculation of ( WNW . The data-point (o) corresponds to the 
numerical result {k^W 1 /it; g) = (3.2; 1.87) of De Raedt et al. 
H TABLE I]. 

In Fig. |] we present a calculation of g as a function 
of k^W /tt for this specific geometry. The result of Do 
Raedt et al. (o) is reproduced but in general the normal- 
ized conductance is seen to have many resonances caused 
by the high aspect ratio of the constriction. In the range 
3 < kpW'/TT < 4, the normalized conductance can be 
anything in the range 1.655 < g < 2 depending on the 
position of the Fermi level and though De Raedt et al. 
]l5[ found the back-focusing phenomena of the Andreev 
reflection to be very robust with respect to changes of 
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the device-parameters, the normalized conductance itself 
certainly depends strongly on the position of the Fermi 
level. The reason is that only those quasiparticles which 
enter the region between the constriction and the inter- 
face can be Andreev reflected and thus contribute to the 
conductance enhancement compared to the normal state 
conductance. 



B. Phase-coherent junction with barrier at interface 

We next consider coherent transport through an NS in- 
terface with a barrier at the interface and a WNW con- 
striction at a distance Li from the interface, see lower 
right insert of Fig. |^. The conductance G NS and the nor- 
mal state conductance G N are found from Eqs. (J3J) and 
(|[) with the transmission matrix in Eq. (|l5|). 
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FIG. 5. Normalized conductance g = Gns/Gn of a coher- 
ent WNW 2DEG-S junction with a barrier as a function of 
the normalized barrier strength Z for k F W' /it = 3.195 (A), 
k F W'/n = 3.2 (o), and k F W'/n = 3.205 (□). The lower 
left insert shows the normalized conductance g as a function 
of k F W' /n for Z = 0. The constriction has an aspect ratio 
Li/W' = 5/1.6, a relative width W/W = 31.72, and the 
cavity has a relative length Lz/W' = 20/1.6. 

In Fig. H we present a calculation of the normalized 
conductance g as a function of the normalized barrier 
strength Z for the device considered by De Raedt et 
al. [y_5|. For the position of the Fermi level (o) consid- 
ered by De Raedt et al., the normalized conductance is 
only weakly suppressed (compared to a system without 
a QPC, see e.g. p9|] ) for low barrier scattering (Z < 1) 
and only for a very high barrier strength (Z > 2) the 
normalized conductance approaches the cross-over from 
an excess conductance (g > 1) to a deficit conductance 
(g < 1). The effect of the barrier for Z < 1 is very sim- 
ilar to the reflectionless tunneling behavior in diffusively 
disordered junctions [^3|,0 where the net result is as if 
tunneling through the barrier is reflectionless. In the case 



of a QPC instead of a diffusive region there is a weak de- 
pendence on the barrier strength and the tunneling is not 
perfectly reflectionless. 

An interesting feature is the non-monotonic behavior 
of g as a function of Z. For Z — > oo, the normalized 
conductance of course vanishes, but in some regions it 
increases with an increasing barrier strength [curve (o)] 
and for Z ~ 1 it has the same value as for Z = 0. This 
is purely an effect of size quantization in the cavity be- 
tween the QPC and the barrier which enters the conduct- 
ing properties because of the fully coherent propagation 
of electrons and holes. However, changing e.g. the po- 
sition of the Fermi level slightly, [curves (A) and (□)], 
changes the quantitative behavior although the overall 
suppression of g with increasing Z is maintained. 



C. Incoherent junction 

In junctions where the propagation in the cavity be- 
tween the QPC and the NS interface is incoherent, the 
so-called sequential tunneling regime, the QPC and the 
NS interface can be considered as two series-connected 
resistive regions pif . This means that 



Gqpc-ns — (G Q p C 



g; 



g~ 
g; 



(31) 
(32) 



where G QPC and G N are found from Eq. (Ji|) with t 



tc 

and t = t 

t = t e . The normalized conductance can be written as 



G 

G c 



QPC-NS 



G N 

~gZ 



G 



QPC 



G N 



G c 



G K 



(33) 



and for W ^ W the major contribution to the resis- 
tance comes from the QPC, i.e. G QPC (G N , G NS ). 
This means that the enhancement of G NS compared to 
G N has a negligible effect on the total conductance so 
that the normalized conductance approaches g ~ 1. 
For an ideal QPC and an ideal interface we have 



G c 



4c' 
h 



M, and G N 



2e- 



M, where 



N is the number of modes in the QPC and M is the 
number of modes at the NS interface. The correspond- 
ing normalized conductance is shown in Fig. ^|. 

The sequential tunneling behavior may provide an ex- 
planation for the unexpectedly small conductance en- 
hancement seen in the experiments of Benistant et al. |T§| ] 
on Ag-Pb interfaces with injection of quasiparticles into 
an Ag crystal through a point contact. The condition for 
the electronic transport to be incoherent is that the dis- 
tance between the point contact and the NS interface is 
longer than the correlation length L c = Min(£j n , Lt), ha 
being the inelastic scattering length and Lt the Thouless 
length j|]6) . For the ballistic device studied by Benistant 
et al. L c = Lt = fiv^/k^T ~ 9/im (at T = 1.2 K) 
which is much shorter than the distance between the 
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point contact and the NS interface (~ 200 /im). Lower- 
ing the temperature will increase the correlation length 
and for sufficiently low temperatures (T ~ 0.05 K) we ex- 
pect a cross-over from the sequential tunneling regime to 
the phase-coherent regime where the Andreev mediated 
conductance enhancement should become observable. 
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FIG. 6. Normalized conductance g = Gqpc-ns/Gqpc-n of 
a QPC NS junction with sequential tunneling through the 
ideal QPC and the ideal NS interface as a function of the ra- 
tio M/N of propagating modes in the QPC and at the NS 
interface. 



scattering length compared to the length of the cavity 

S 

For the sequential tunneling regime we find that the 
conductance enhancement vanishes as the number of 
modes at the interface becomes much larger than the 
number of modes in the QPC. 

Our calculations show that the S-matrix approach 
provides a powerful alternative to the time-dependent 
Bogoliubov-de Gennes approach of De Raedt et al. jl5| in 
detailed studies of the conducting properties of nanoscale 
2DEG-S devices. Even though the back-focusing phe- 
nomenon of the Andreev reflection is robust against 
changes in the geometry p5| , the electron- hole conver- 
sion efficiency itself is not. 

Finally, we stress that for a quantitative comparison to 
experimental systems, it is crucial to take different Fermi 
wave vectors and effective masses of the 2DEG and the 
superconductor into account pfjfl. 
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VI. DISCUSSION AND CONCLUSION 



For an ideal 2DEG-S interface with a QPC in the nor- 
mal region, the normalized conductance g = G NS / G N de- 
pends strongly on the position of the Fermi level and 
only when the Fermi level corresponds to a conductance 
plateau a doubling of the conductance is found. The de- 
viations from the factor-of-two enhancement, when the 
Fermi level does not correspond to a plateau, can be sig- 
nificant and for a particular example of the WNW con- 
striction we find that the normalized conductance can be 
suppressed to g ~ 1.5 in a system with only two prop- 
agating modes in the constriction. In the presence of a 
barrier at the 2DEG-S interface, the normalized conduc- 
tance depends strongly on the longitudinal quantization 
in the cavity set up by the QPC and the barrier. Depend- 
ing on the barrier strength, the length of the cavity and 
the position of the Fermi level, this longitudinal quanti- 
zation may give rise to both constructive and destructive 
inferences in the transmission and thus also in the con- 
ductance. Perhaps surprisingly, the effect of the barrier 
is very much suppressed (compared to a system without 
a QPC, see e.g. due to a very strong back-scattering 
at the return of the quasiparticles to the normal probe. 
The localization of quasiparticles in the cavity gives rise 
to an almost reflectionless tunneling through the barrier 
as it is also found in systems with a diffusive normal re- 
gion JT^JT^]. The interferences due to localization in the 
cavity will be smeared by a finite temperature and they 
are also expected to be suppressed by a finite inelastic 
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